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ABSTRACT 



Context. For a long time the consensus has been that star formation rates are higher in the interior of spiral arms in 
galaxies, compared to inter-arm regions. However, recent studies have found that the star formation inside the arms 
is not more efficient than elsewhere in the galaxy. Previous studies have based their conclusion mainly on integrated 
light. We use resolved stellar populations to investigate the star formation rates throughout the nearby spiral galaxy 
NGC 5236. 

Aims. We aim to investigate how the star formation rate varies in the spiral arms compared to the inter-arm regions, 
using optical space-based observations of NGC 5236. 

Methods. Using ground-based Ha images we traced regions of recent star formation, and reconstructed the arms of the 
galaxy. Using HST/ACS images we estimate star formation histories by means of the synthetic CMD method. 
Results. Arms based on Ha images showed to follow the regions where stellar crowding is higher. Star formation rates 
for individual arms over the fields covered were estimated between 10 to 100 Myr, where the stellar photometry is 
less affected by incompleteness. Comparison between arms and inter-arm surface star formation rate densities (T>sfr) 
suggested higher values in the arms (~0.6 dex). Over a small fraction of one arm we checked how the ^sfr changes for 
the trailing and leading part. The leading part of the arm showed to have a higher T>sfr in the age range 10-100 Myr. 
Conclusions. Predictions from the density wave theory of a rapid increase in the star formation at the edge where the 
stars and the gas enter the density wave are confirmed. The T>sfr presents a steep decrease with distance from the 
center of the arms through the inter-arm regions. 

Key words, galaxies: Individual - NGC 5236 galaxies: Star formation 



1. Introduction 



Studies of OB stars (see e.g. [Morgan et al.lll953t iMcGruderl 
119751 : iMuzziol [l97l iKaltchev a 2009), show that the con- 
centration of these stars is higher in the Sagittarius arm 
of the Galaxy. This suggests an active star formation in 
the arms of spiral galaxies. The creation of spiral arms is 
a problem that has been studi ed since the late 60's (e.g. 
iLin fe Shulfl965 iRobe rts 1969). The density wave theory 
explains how the spiral arms can be formed and can re- 
main stable over time for an isolated galaxy. The theory 
predicts more active star formation in the arms, where 
the gas compression induced by the density waves triggers 
the process (see e.g. [Lin fc Shul 11964 iBash fc Vissedll98lt 
Kna pen et atlll992l 119961 : iKurtz et al.ll2002t iGrosbol et all 



20061 ). Alternative theories to explain the spiral arms in 



galaxies have been proposed. In the modal theory, inward- 
moving waves reflect or refract off at the center of a galaxy 
(even in the galaxies with a bar in the center), and then 
the w ave comes back out as leading or trailing spira l arms 
(e.g. iMarkl ll974 _ lLau fc Markl Il976t iBertin et all 119891 : 
El megreen et al.lll992l for theoretical and observational ap- 
proaches). Another theory, know as the Stochasti c Self- 
Propagating Star Formation ([Mueller fc Arnettll 19761 ), sug- 
gests that episodes like supernovae, shock wave or gravita- 
tional interactions are responsible to propagate and trigger 
star formation (e.g. lGerola fc Seiden 1978; Seid en fc Gerolal 



ll9T9b iFeitzinger et al.l Il98ll for theoretical and observa- 
tional approaches). This theory is very good at explain- 
ing flocculent galaxies, while not grand-design ones. It is 
possible that the whole process is a combination of these 
theories. However, in this paper, we will test only results 
expected from the density wave theory. 

If the density wave theory is correct, the process of 
gas being compressed by the waves should lead to many 
observable effects, e.g. different star formation rates or 
color-gradients across the arms. Assuming a constant an- 
gular velocity of the spiral density wave and an approx- 
imately flat rotation curve of the stellar component, in- 
side the corotation radius the gas will overtake the density 
wave, which will produce an increase in the star forma- 
tion when compressed. Outside the corotation radius the 
wave overtakes the gas. Consequently, the stars will drift 
and age, creating a color gr adient that can be observed. 
Ma rt mez- Garcia et al.l ([2009) have studied the color gradi- 
ent across spiral arms of 13 spiral galaxies. In their work, 
a reddening free index was used to study this process, con- 
cluding that azimuthal color gradients are common in spiral 
arms of disk galaxies. Previous works have tri ed to relate 
the sp i ral arms with the star fo r mation (see e.g. Allen et "aD 
1986; Tilanus fc Ailed Il993t IStedman fc Knapenl 12001 



Knapen et al.l 120101 : ISanchez-Gil et al.l 120111 ). These stud- 



ies have shown the relation between the star formation and 
the spiral structure using different gas components over dif- 
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ferent parts of a large set of galaxies (e.g. M83, M51, M100, 
M101, among others). 

Compression of the gas across the galaxy is expected to 
be observed as regions with enhanced star formation. The 
(common) components used to estimate the star formation 
efficiencies (SFE) in a galaxy are the gas and the stellar 
population, the latter being commonly measured through 
their integrated light. When H 2 i s use to trace star form- 
ing regions (e.g. iFovle et al.ll2010[ h there seemed to be no 
specific correlation between the SFE and the spiral arms. 
However, studies that estimate the SFE as the fraction of 
Ha to CO and/or HI maps do find a correlation between 
the SFE and the spiral arm (e.g. lLord fc Youngl Q.990; 
ICepa fc Beckmanlll990l : iKnapen et al.lll992l ). Nevertheless, 
it is important to note that most of the studies of star for- 
mation (either efficiency or rate) and their relation with 
the arms in spiral galaxies have been done using gas com- 
ponents, namely HI, H2, CO, etc combined with SFR trac- 
ers (e.g. Ha), and unresolved stellar pop ulations, mea- 
sured through their in tegrated light (e.g. IKnapen et al.l 
1996J iLerov et al.ll2008l ). Tracers based on gas components 
are assumed to indicate the regions where SFE is high, 
which could lead to erroneous estimations, as stated by 
IFovle et al.l (|2Q1Q[ ). 

As a new approach, we present in this paper the use of 
resolved stellar populations as a tool to study how the sur- 
face star formation rate density (Tisfr) varies in the arms 
and inter- arm areas of the nearby spiral galaxy NGC 5236. 
This galaxy has been selected because it is one of the clos- 
est face-on, nearby, grand-design spiral galaxies, where spa- 
tial resolution allows a study of this kind. The low inclina- 
tion in the line-of-sight reduces the effect of internal red- 
dening. Previous estimations of t he T^sfr show values of 



13 x lQ- 3 M^yr- 1 Kpc~ 2 (e.g. [Larsen fc Richtlerl |2QQQ|; 
ICalzetti et all l2Q10t ISilva- Villa fc Larsenl [2QlTh . which is 
high compared to oth er normal spiral galaxies (e.g. 
lLarsen fc Richtlerl I2QQQL see table 2). A large number of 
studies have been done in the center of this galaxy, study- 
ing different properties through different wavelengths, all 
of them suggesting an increasing activ ity in the star for- 
mation during the last ~ lQMyr (e. g. Rvder et al. 1995; 
i Harris et al1l200lUHoughton fc Thatt e 2008; Knapen et all 



120101 ). Activity in the center of the galaxy combined with 
the high levels of T<sfr in the disk, show that NGC 5236 is 
still actively forming stars, making it ideal to study differ- 
ences in the star formation and the possible relations with 
environment and location. 

Kinematic s t udies of NGC 5236 were presented by 
iLundgren et al.l (|2004l ) based on gas (H2+HI) measure- 
ments. The estimation of the gas surface density in the 
arms done by Lundgren et al. is higher than the Toomre's 
value for stability (Q oc S" 1 ^). This is possibly causing in- 
stabilities in the arms of the galaxy, potentially leading to 
star formation. However, their estimation of the gas sur- 
face density in the inter-arm r egions do not show the same 
high values. In a further work. [Lundgren et al.l (|2008l ) used 
far-UV, B and Ha integrated light to estimate star forma- 
tion rates, while CO was used for gas maps. Lundgren et 
al. conclude that the star formation presents higher levels 
in the nuclear regions, close where the bar ends, and in the 
arms of the galaxy. The authors also found an increased 
SFE along the arms of this galaxy. 

This paper is structured as follows. We introduce the 
observations, optical and Ha, used to study the field stellar 



population and the regions of recent star formation, respec- 
tively, in Sect. 2. Section 3 is devoted to review how the 
optical bands were used to run the photometry of the field 
stars. Using Ha (i.e. tracing recent star formation), we will 
introduce in Sect. 4 the method used to re-construct the 
arms of the galaxy. Using the photometry from Sect. 3, we 
estimated the star formation history of different groups of 
field stars in Sect. 5. Finally, we present our discussion and 
conclusions in Sects. 5 and 6, respectively. 

2. Observations 

We used images from Hubble Space Telescope and Cerro 
Tololo observatory, covering different optical wavelengths 
and Ha. 

2.1. BVI observations and data reduction 

We use observations of NGC 5236 taken by the Hubble 
Space Telescope (HST), using the Advanced Camera for 
Surveys (ACS). The instrument has a resolution of 0'/05 
per pixel. With a distance modulus of 28.27 (~ 4Mpc, 
iThim et aDl2QQ3h . 1 pixel corresponds to ~1 pc in our im- 
ages. 

The images of this galaxy have been taken in the optical 
bands F435W (~ £), F555W (~ V), and F814W (~ I), 
with exposure times of 680 sec for the bands B and V, and 
430 sec for the band I. Our observations were taken in 2004 
as part of Cycle 12, centered at a : 13 : 37 : 00 and S : — 29 : 
49 : 38 and a : 13 : 37 : 06 and 5 : -29 : 55 : 28 (J2000) 
for the first and second field observed, respectively. Figure 
[T] presents a DSS image indicating the regions observed. 

The standard STScI pipeline was used for the initial 
data processi ng. ACS images were d rizzled using the mul- 
tidrizzle task ([Koekemoer et al"1l2002f ) in the STSDAS pack- 
age in IRAF using the default parameters, but disabling 
the automatic sky subtraction. Object detection for field 
stars was performed on an average B, V, and I image, us- 
ing daofind in IRAF. 

2.2. Ha observations 

We use an archival Cerro Tololo Inter- American 
Observatory (CTIO) 1.5m telescope image of NGC 5236, 
taken in 2006 as part of the Survey for I onization in 
Neutral Gas Galaxies fSINGG. IMeurer et all I2QQ61P L The 
survey used the Ha and R bands, with a resolution of 
0'/4 per pixel. The total exposure time of the observation 
for NGC 5236 was 1800 sec, with center coordinates at 
a : 13 : 37 : 02 and S : -29 : 52 : 06 (J2000). The image 
used in this paper is based only on the Ha filter. For 
det ails on the observa tions and calibration of the image 
see IMeurer et all (|2006[ ). 



3. Field stars photometry 

Details of our method o f analysis can be found in 
ISilva- Villa fc Larsenl (|201Q[ ). Below we will reiterate the 
main points of the procedures we used to carry out pho- 
tometry on our data. 



1 SINGG is a subsample o f the HI Parkes All Sky Survey 
fHIPASSl. lMever et all (|200l 
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Fig. 1. M83 DSS image. Blue lines delineate the ACS point- 
ings used in this paper. 
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Fig. 2. Left column: Color-magnitude diagrams for the Arm 
1, Arm 2, and inter- arm regions. Right column: Histogram 
of colors of main sequence stars inside the boxes marked 
over the CMDs. Errors are Poissonian. Vertical dashed lines 
represent the peak of the Gaussian distributions. 
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Fig. 3. Upper panel: Blurred image used for the detection of 
pronounced Ha regions. Lower panel: Original CTIO image 
of NGC 5236. Overplotted in each panel are the estimated 
locations of the arms. In the lower panel the two fields ob- 
served with the HST/ACS are marked as yellow regions. 
Gree n dot- dashed line is th e corotation radius located at 
170" (jLundgren et al.ll2004ft . Dash-dotted lines in blue de- 
note annuli at 200 / /0, 140'/0 and llO'/O, see text for details. 



Due to the crowding, we performed PSF photometry for 
field stars. Using a set of bona-fide stars visually selected 
in our images, measuring their FWHM with imexam, we 
construct our point-spread function (PSF) using the PSF 
task in DAOPHOT. This procedure is employed in the same 
manner for each band (i.e. B, V, and I). The PSF stars are 
selected individually in each band, in order to appear bright 



and isolated. PSF photometry is done with DAOPHOT in 
IRAF. 

Our PSF-fitting magnitudes are corrected to a nominal 
aperture radius of 0'/5, following standard procedures. From 
this nominal value t o infinity, we apply the corrections in 
ISirianni et all (|2QQ5f ). 
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Fig. 4. Surface density of stars for the ACS fields. Red lines 
represent the path of the arms. 

HST zero-pointfl were applied to the PSF magnitudes 
after applying aperture corrections. The zero-points used 
in this work are ZP B = 25.77, ZP V = 25.72 and ZPj = 
25.52 magnitudes. Typical errors of o ur photometry do not 
chang e dramatically from the ones in [Silva- Villa fc Larsenl 
(|2Q1q[ see its Fig. 2). 

The final color-magnitude diagrams (CMD) for the stars 
in the arms and the inter- arm regions (see Sect. 5 for defi- 
nition of the stars that belong to the arms and in the inter- 
arm region) are presented in Fig.[2j left column. In the same 
figure we investigate whether there are any indications of 
differences in the mean extinction from one region to an- 
other. We created histograms of the main sequence stars, 
defined to be stars in the color range —0.3 < (V — I) < 0.7 
and in the magnitude range — 10 < My < —5, as indicated 
by the red boxes over the CMDs (see Fig. [2j right column). 
The histograms were made with variable bin widths using 
100 stars per bin. The number of stars is normalized to the 
width of the bin. We fitted Gaussians to the observed distri- 
butions to estimate the maximum of the distributions. We 
do not observe a large variation among the Arm 2 and the 
inter- arm areas, where the peak of the distribution is close 
the same value (~ 0.10±0.02). However, the Arm 1 presents 
a slight shift in the peak of the distribution (~0.13±0.02) 
compared to Arm 2 and inter- arm regions. This difference 
is not large (< 0.03, and errors overlap), but we note that 
the Arm 1 is closer to the galactic center, where extinction 
can be affecting the observations. There is also an apparent 
increase in the distributions close to (V — I) ~0.5. Inside 
our photometric errors, separating the main sequence stars 
from the blue He burning phases is not straightforward, but 
the count of stars can give an indication of the presence of 
these stars, as seen in Fig. [2j where at (V — I) « 0.5 the 
distribution of stars clearly deviates from the Gaussian fit. 

4. Defining the spiral arms over NGC 5236 

To identify t he arms of NGC 5236, we followed the method 
described bv lScheepmaker et all ([20091 ), using Ha as indica- 
tor. The CTIO image of the galaxy was cropped to remove 
as much background as possible. Over the new (sub) image, 
we use a Gaussian kernel (with a 20 pixels sigma) to blur 
the image, enhancing the regions where Ha appears to be 
more concentrated. Because of the high concentration of 
gas in the center of the galaxy, we manually mask this part 

2 www.stsci.edu /hst / acs / analysis / zeropoints / #tablestart 



of the image, which will allow better analysis of the re- 
gions outside of the center. We analyzed the blurred image 
with Daofind in IRAF to find the places where Ha is more 
concentrated (minimum data value of 400 counts over the 
background). The coordinates retrieved by IRAF were vi- 
sually inspected to remove unnecessary detections, i.e. de- 
tections that could affect the estimations of the arm's path 
and that are not part of the ACS field-of-view. We used 
a cubic spline interpolation to fit the arms of NGC 5236. 
Figure [3] depicts the CTIO images. The upper panel is the 
cropped and blurred image, overplotted with the estimated 
location of the arms. The lower panel shows the original 
CTIO image, overplotted with the location of the arms (red 
lines) and the two fields covered by the ACS observations 
(yellow squares). Conversion of the coordinates from the 
CTIO coordinate system to the ACS coordinate system was 
done using wcsctran in IRAF. First, using the header of the 
CTIO image, we converted the CTIO coordinate system to 
WCS coordinates, and then, using the same procedure, we 
moved finally to the ACS coordinate system. The different 
circumferences in Fig. [3] are marking the distances 200'ZO, 
140'/0, 170'/0 and llO'/O used to estimate the velocity of a 
particle at different radii (see Sect. 6). 

We create the density plot of the resolved stellar pop- 
ulations over the ACS fields with the arm's path cre- 
ated using the Ha. After overplotting the arms, we ob- 
serve that the arms follow the region where the density 
of stars is higher, as seen in Fig. SJ From the same fig- 
ure, it is observed that other regions have high density 
of stars, e.g. in field 2 there is a possible "feather" be- 
tween the two arms, and in field 1 we observe many re- 
gions with high densities of stars, which is expected, due 
to the proximity to the center o f the galaxy, where the 
concentration of gas is higher (e.g. ICrosthwaite et al.l l2002; 
iLundgren et al.l l2004f ). There is large observational evi- 
dence related with "feathers" and/or "spurs" close (or at- 
tached) to spir al arms, which is supported by theoretical 
studies (see e.g. IScoville et aH 120011: iKim fc O striker 2002; 
IShettv fc Ostriker]l2006l : lLa Vigne et al.ll2006f ). " 

For the rest of this paper we will refer as Arm 1 to 
the arm that is fully inside the corotation radius (see green 
dotted line in Fig. [3]), while the Arm 2 has one part inside 
the corotation radius and one part outside of it. Figure [3] 
pr esents the two arms and the corotation radius, estimated 
bv ILundgren et all ([2004D to be located at 170". 



5. Selection of stars and star formation histories 

We aim to study here how the star formation varies from 
the center of the spiral arms to the inter-arm regions. We 
use the resolved stellar population detected over our ACS 
fields and estimate their star formation histories at a fixed 
distance from the center of the arm. 

5.1. The arms 

Stars were selected every 0.2 Kpc, assuming that the arms 
presented above mark the "center" of the distribution. 
Figure [5] shows the distribution of the stars which follow 
the arms at the distances 0-0.2, 0.2-0.4, and 0.4-0.6 Kpc, 
represented with the colors yellow, green, and blue, respec- 
tively. In our second field the distance between the two 
arms would allow to cover a larger area, however, we kept 
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Fig. 5. Selection of the resolved stellar population over the arms of NGC 5236. Red lines show the arms of the galaxy 
based on Ha images. Yellow, green, and blue represent the selection of stars at different distances that belong to the 
arms. Orange regions represent the inter- arms. 



the same distances as in the first field in order to have com- 
parable results among the fields. 

For each of the selected group of stars, star forma- 
tion histories ( SFH) were calc ulated using the synthetic 
CMD method ([Tosi et al.lll99lh . A description and tests 
of the IDL-based program used t o estimate the SFH is 
given m ISilva- Villa fc Larsenl (1201 Oh . For more applications, 
see ISilva- Villa fc Larsenl (l201lt). W e have updated the code 
presented by ISilva- Villa fc Larsenl (|2010l ) allowing the use 
of multiple metallicities. When using multiple metallicities, 
the code models the Hess diagram as a weighted sum of all 
isochrones in the input library. The final SFH is obtained 
adding together the individual SFR over time for each as- 
sumed metallicity. 

The parameters used for the estimation of the SFH in 
NGC 5236 are: Distance modulus of 28.27, solar and LMC- 
like metallicities, a binary fraction of 0.5 with a mass ra- 
tio between 0.1-0.9 (assuming a flat distribution and no 
binary evolution), and using the color combination V-I. 
We normalized our estimations by the areas covered, hav- 
ing then the surface star formation rate density (T,sfr 
[M0yr -1 Kpc -2 ]). The estimation of the different areas cov- 
ered was carried out following a similar procedure as for the 
selection of the stars. We create an image of the size of the 
ACS fields and calculate the distance of each pixel to the 
arms. Having the total amount of pixels under the desired 
area, the areas were calculated following the relation: 



Motal X 



Npix Al 
Npix Atotal ' 



(1) 



where A^ is the area to be calculated, A tota i is the total area 
covered by the ACS (21.3 Kpc 2 , each field), and Npix^. and 



Npix^ totaZ are the number of pixels for the area i and the 
total number of pixels, respectively. 

Figure [6] illustrates the fits done for field 1 for the three 
areas covering Arm 2. We note that the fits are far from 
perfect. Trying to understand what could be the cause of 
the mismatch between observed and modelled CMDs, we 
repeated the fits varying the input parameters. As an exam- 
ple, the right column in Fig.[6]shows the fits using only solar 
metallicity (instead of solar and LMC). It is clear that solar 
metallicity leaves a "gap" between the main sequence stars 
and the blue loop, which is not seen in our observed CMDs 
(first column). The combination of both metallicities gives a 
better fit to the data, as the blue loops extend to higher ef- 
fective temperatures (bluer colours) for LMC-like metallic- 
ity, and we can not reject the possibility that higher metal- 
licities are present. Regardless of the improvement when 
using a combination of metallicities, we note that neither 
red nor blue He burning phases are well fitted by our pro- 
gram, appearing bluer and redder, respectively, in our final 
fits, when compared with the observations. Assumptions of 
no binaries did not show any improvement. Following the 
small change in extinction, as suggested in Sect. 2, we used 
a time dependent extinction, but we did not see any im- 
provement. The middle column presents the best fit which 
combines the parameters described in the previous para- 
graph. 

In view of the difficulties reproducing the observed 
CMDs, caution should clearly be exercised when interpret- 
ing SFHs inferred from these fits. 

For each field and each arm, between 10 and 100 Myr 
(age range less affected by incompleteness), we estimated 
the mean values for ^sfr- We also combined results for 
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Fig. 6. Left colum: CMDs for the Arm 2 in the field 1 at dif- 
ferent distances from the center of the arm. Middle colum: 
Best fit CMDs assuming solar and LMC like metallicities, 
a binary fraction of 0.5 with a mass ratio between 0.1-0.9, 
and a fixed A# foreground extinction of 0.29. Right colum: 
Best fit CMDs assuming solar metallicity, a binary fraction 
of 0.5 with a mass ratio between 0.1-0.9, and a fixed A# 
foreground extinction of 0.29. The fits were done using the 
red boxes shown in the figure. 



each arm, adding the results from each field, see Fig. [7] and 
[8j We summarize our results in table [TJ 

Independent measurements of T<sfr (i-e. per field, per 
arm, per area covered) are shown in Fig. [71 where Arm 
2 suggest that during the past 100 Myr there has been an 
increase in the T^sfr (right panels). As the individual SFHs 
are quite noisy, we joined the results for all areas of Arm 2, 
but for the Arm 1 we only used the results from the field 2, 
because the areas covered in the field 1 were too small and a 
clear separation between the arm and the inter-arm region 
was not easy to performed, compromising the results. The 
results are presented in Fig. [51 From this figure it becomes 
clear that there is an increase in the T<sfr between 10 to 
100 Myr for Arm 2. This increase is observed for each one 
of the areas covered in that particular arm, which suggest 
a physical effect. In Sect. 6 we interpret this result. 

5.2. The inter-arm regions 

To estimate the T<sfr for the inter- arm regions we used 
the stars outside the arms (orange regions in Fig. [5]). Stars 
selected in this way were passed to our program and SFHs 
for each field were estimated (see Fig. [8] right panel). The 
areas were calculated as the total area (A tota i) minus the 



^sfr and Arms over NGC 5236 

Table 1. Surface star formation rate densities at different 
distances from the center of the arms. Errors are the stan- 
dard deviation of the mean. Arm 1 values are for field 2 
only. 



FiekLArm 


Dist arr n 


Area 


^SFR 




[Kpc] 


[Kpc 2 ] [xl(T 


^Moyr^Kpc- 2 ] 






Arms covered 




Arml 


0.0-0.2 


1.06 


67±24 


Arml 


0.2-0.4 


1.05 


43±14 


Arml 


0.4-0.6 


1.01 


29±18 


Arm2 


0.0-0.2 


3.29 


67±15 


Arm2 


0.2-0.4 


3.44 


41±19 


Arm2 


0.4-0.6 


3.54 


31±17 




Inter- arms covered 




Fieldl 




12.84 


23±11 


Field2 




15.24 


8±3 



areas covered by the arms, then having an estimation of 
the T>sfr (see table [I]). We added the independent esti- 
mations for the two fields and plotted the mean value in 
Fig. [9] as a continuum value, which is lower than the values 
derived for the arms, making a clear distinction between 
the SFH's in the arms compared to the inter-arm areas. 
It is important to note the differences among the fields. 
The first field, which is close to the center of the galaxy 
displays a higher T<sfr, even suggesting a small increase 
over time. The second field, which is located further away 
from the center, displays a lower T>sfr- This separation in- 
dicates differences that could be related with the amount 
of material that could potentially be used to form stars. It 
is expected to find more gas concentrat ed in the center of 
the galaxy than in the oute r parts (e.g. ICrosthwaite et all 
20021 [Lundgren et al.ll2004l ). Also important to note is the 
increase between ~20-100 Myr. This could be related to the 
limits adopted for the arms and inter-arm regions, where 
the increase in the number of stars at the " edge" of the arms 
could be contaminating the inter- arm areas. Also, as noted 
in the previous sections, there are regions in the inter-arms 
that present high levels of Ha, suggesting recent star forma- 
tion activity. Nevertheless, keeping in mind the discussion 
in the previous section, the increase could be an artifact of 
the fitting process. 

6. Discussion 

The Tisfr of the arms in NGC 5236 estimated here suggest 
that the arms of the galaxy have increased their star forma- 
tion rate during the last 100 Myr by ^0.2 dex (see Fig. [8j 
left and middle panels). In the right panel of the same fig- 
ure we present the T>sfr of the inter- arms, which presents 
three important features: (1.) There is a clear difference in 
the star formation history among the two fields, where the 
field 1 presents a higher T>sfr in comparison with the field 
2 (see table [IJ; (2.) there is an apparent increase in the star 
formation history, which could be due to different factors, 
e.g. feather/spurs with active star formation or small differ- 
ences in the separation between arm and inter-arm areas; 
and (3.) both fields present a mean T>sfr lower than in the 
arms areas (see table H}. 

Figure [9] shows the variation of the mean T>sfr (be- 
tween 10 and 100 Myr) as a function of distance from the 
center of the arms, and the T*sfr for the inter- arm re- 
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gions, as defined in the previous section. Starting from the 
innermost region of the arms, traced by Ha, we observed 
a very fast and steep decrease moving outward through 
the arms. The estimations of T^sfr show a clear difference 
between the arms and the inter-arm regions, the former 
having a higher value in the center (by ~ 0.6 dex), while 
reaching similar levels at the external regions, where the 
separation between arm and inter- arm becomes dificult. In 
Fig. [4] we showed that Ha traces the regions where the 
crowding of stars is higher, suggesting that the recent star 
formation is higher in the arms of the galaxy. There are 
other regions in the same figure where large concentrations 
of stars are present, but based on the procedures described 
here, in those regions Ha could be interpreted as "feath- 
ers" and/or "spurs" based on our procedures. We do not 
disregard the possibility of having recent star formation in 
the inter- arm regions, pos s ibly in duced by e.g. supernovae 
explosion. iLundgren et al.l (|2QQ8f ) suggested that in the bar 
and in the arms of M83, the T*sfr and the SFEs are higher 
than in the inter-arm areas. Our measurements of T<sfr 
showed to be higher in the arms of the galaxy, in agree- 
ment with the work done by ILundgren et all (|2008f h when 
compared to the inter-arm regions (at least qualitatively), 
however, we do not draw any conclusion regarding the bar 
of the galaxy. 

We tested results expected from the density wave theory 
studying the T,$fr in the fields covered in this paper. Using 
increasing distances from the center of the arms (as defined 
in this paper through Ha) a comparison was made to see 
the spatial variations in the arms of the galaxy, looking for 
the differences between the arms and the in t er- arm regions. 

The density wave theory (|Lin fc Shul Il964t iRobertsI 

Il969h explains how it is possible to have semi-stable arms 
in spiral galaxies. We test predictions from this theory us- 
ing resolved stellar population. If the gas is overtaken by 
the density wave (or viceversa, depending on the radius), 
then, it is expected to observe higher star formation rate per 
unit area at the edge where the gas enter s the density wave 
due to the compression of the gas (e.g. I Allen et all 1 19851 : 



Elmegreen 2011). Also, from the density wave theory it is 
expected to observe young stars close to the regions where 
the gas is more concentrated (i.e. where the star formation 
is happening), but with time and rotation, stars evolve and 
travel through the arm. This process can be observed as 
a color-gradient across the arms. Theoretical predictions 
of the time ranges for color -gradients were estimated by 
Ma rtmez-Garcfa et al.l ([2009) using the reddening-free pa- 
rameter Q(rJgi) (where r, J, g, i are the photometric bands 
used). Martinez- Garcia et al. show that the maximum of Q 
is reache d at ~ 2 x 10 7 yr after a burst of stars (assuming a 
Sal peterl (|l955f ) IMF), and then starts to decline, reaching 
lower values at ~ lOOMyr after the burst. This estimate 
is of course affected b y inclination, the photometric bands 
used, metallicity, etc. Mar tmez-Garcfa et al.l ([2009) stud- 
ied the color gradient variation over different spiral galax- 
ies. The authors used a sample of 13 spiral galaxies, and 
showed how a fraction of the arms analyzed clearly showed 
the expected color gradient suggested by the density wave 
theory. We used resolved stellar populations over the Arm 
2, in field 1, where the statistics are better (see Fig. [5]), 
to observe the variations of the T>sfr inside this (part of 
the) arm. Separating the areas at 0.2-0.4 Kpc and 0.4-0.6 
Kpc into "trailing" and "leading" parts of the arm (left 
and right in Fig. |5j), we est imated T>sfr separate ly. Using 
the rotation curve found in ILundgren et ahl (|2004l ). we esti- 
mated how much time a particle will need to move 0.2 Kpc 
at different radii. Ins ide corotation radius (Rcr = 170", 
ILundgren et al.| [2004. ~3 Kpc for the distance modulus as- 
sumed in this paper), the time for a particle to move that 
distance is close to ~5 Myr. Outside corotation, the time is 
increased by a factor of ~3 times. As observed in Fig. [3] the 
circumference at different radii are not crossing perpendic- 
ular to the arms. Because of this, the real path a particle 
will cross is larger than 0.2 Kpc. We decided to estimate 
mean T>sfr every lOMyr, and check how it varies with time. 
Figure [10] shows the T>sfr as a function of position in the 
arm for different time intervals. Inside corotation, the gas 
overtakes the density wave, creating a very steep increase in 
the first few Myr, which must start to decrease with time, 
as observed in the figure. Fig. [TT1 shows the average of the 
^sfr over an age range between 10-100Myr. As indicated 
by the arrow, and being (mostly) inside corotation, a very 
steep increase in the T<sfr is observed (by a factor of ~0.4 
dex), reaching a maximum where the location of the Ha is 
denoting the center of the arm, and then showing a decrease 
with time. The small fraction of the arm outside corotation 
was not analyzed separately due to the low statistics and 
the lack of smoothness of the arm at that location (see Fig. 

13). 

As a consistency check to our results, we selected the 
massive bright stars over the CMD for Arm 2 in field 
1, defined as stars between —0.3 < (V — I) < 0.7 and 
— 10 < My < — 7. The number of stars in each area were 
divided by the area covered and then passed to logarithmic 
scales. The estimated stellar densities were shifted to the 
same levels of the SFH estimations for the comparison (see 
Fig. [10] red dash-dotted line). We observed that the dis- 
tribution of these stars is in excellent agreement with the 
SFRs derived from the fits to the CMDs. 

Thes e results are in good agreeme nt with the conclu- 
sions of Martf nez-Garcfa et al.l ([20091 ) for spiral galaxies, 
and with the predictions of the density wave theory, which 
states that the leading part of the arm should have a higher 
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star formation compare to the trailing part of an arm, inside 
the corotation radius. 



7. Summary and conclusions 

Using Ha CTIO images of the galaxy NGC 5236, we 
have traced the arms of the galaxy. Superb resolution of 
HST/ACS observations have allowed us to identify the stel- 
lar populations, which we have used to observe the varia- 
tions of T<sfr in the inner parts of the arms and compare 
how this parameter behaves with increasing distance from 
the center of the arms through the inter- arm regions. 

Ka showed to be a good tracer of higher concentrations 
of stars across the regions covered with our HST obser- 
vations. Star formation histories were es timated using the 
synthetic CMD method ([Tosi et al .11 199 it ) across the arms 
and the inter-arm regions separately. However, the fits to 
the CMDs done are far from perfect. Verifying different 
parameters did not further improved our results, however 
we observed that the inclusion of more than one metallicity 
helps to reproduce better the observed CMDs. We observed 
that the arms of NGC 5236 present a higher T<sfr, com- 
pared to the mean value for the inter-arm regions combined 
by ^0.6 dex. Our analysis of r esolved stellar populat ions 
leads to similar conclusions as iLundgren et al.l ([20081 ) for 
this galaxy, which were based on gas, dust, and integrated 
light studies. 

Comparison of the trailing and leading part of the Arm 
2 (only for the observations in field 1) indicates that the 
region leading has a higher T,sfr. Th i s resu lt is in good 
agreement with iMartmez-Garcfa et al.l (|2009f ) , who used a 
reddening- free photometric index ( Q) to trace the variation 



Fig. 11. Variation of T^sfr across the Arm 2, observed in 
field 1. Variation averaged over the whole age range 10-100 
Myr. Errors are the standard deviation of the mean. 



of the stars formation across a spiral arm, using a sample 
of 13 spiral galaxies. Both results are in agreement with the 
density wave theory, which suggested that the region hit- 
ting, either the gas or the spiral density wave, must induce 
a higher star formation due to the increase in the density 
of the gas. 

An alternative to the spiral density wave theory is 
the theory of stochastic self-propagating star formation, 
whi ch is quite successful a t explaining flocculent galax- 
ies ([Mueller fc Arn ett 1976). In this theory, star formation 
propagates through shock waves produced by supernova 
events. These shock waves induce new regions of active star 
formation, which combined with the differential rotation of 
the galaxy, create the fragmentary spiral patterns observed 
in flocculent galaxies. However, it is unlikely that this sce- 
nario can account for grand-design spiral galaxies such as 
M83. It is possible that supernova shocks and other feed- 
back leads to some propagation of star formation within 
the spiral arms, but our observations are not well suited for 
addressing this issue as we cannot easily age date individual 
stars at the very young ages involved. Also, our broad-band 
imaging does not allow us to observe the detailed structure 
of the gas, which could be a good indicator of this process. 
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